Anisotropy and Penetration Depth of MgB 2 from n B NMR 
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The n B NMR spectra in poly crystalline MgE>2 were measured for several magnetic fields (1.97 T 
and 3.15 T) as a function of temperature from 5 K to 40 K. The composite spectra in the super- 
conducting state can be understood in terms of anisotropy of the upper critical field, jh, which is 
determined to be 5.4 at low temperature. Using Brandt's algorithm 1 the full spectrum, including 
satellites, was simulated for the temperature 8 K and a magnetic field of 1.97 T. The penetration 
depth A was determined to be 1, 152 ± 50 A, and the anisotropy of the penetration depth, 7a, was 
estimated to be close to one at low temperature. Therefore, our findings establish that there are 
two different anisotropies for upper critical field and penetration depth at low temperatures. 



The discovery of unusually high superconductive tran- 
sition temperatures of MgB2, a simple bimetallic com- 
pound superconductor—, has attracted considerable in- 
terest from both theory and experiment. Reconsidera- 
tion and extension of BCS theory to two-band super- 
conductivity has successfully accounted for experimental 
observations 3 ^! 5 -^^. Nonetheless, the relation between 
anisotropy of the upper critical field and the penetration 
depth is still a controversial issue. Generally there are 
two points of view. One holds that there exist two dif- 
ferent anisotropies at low temperatures, and 7a, for 
upper critical field and penetration depth respectively. 
They have different temperature dependence and merge 
at a common value at T,£&1£. The other perspective 
is that there is only one anisotropy parameter, and it 
is field dependent* 5 - ' 11 ! 12 . Moreover, small angle neutron 
scattering (SANS) gives different results on the penetra- 
tion depth anisotropy on single crystal and powder MgB2 
samples ^ 13 ' 14 . In general, it is more of a challenge to de- 
termine the absolute value of the penetration depth as 
compared with its temperature dependence. Although 
muon spin resonance (/zSR)^ 5 - and nuclear magnetic res- 
onance (NMR) methods 1 ^ have often been used to obtain 
an absolute value of the penetration depth, the applica- 
tion of these resonance techniques to determine the pen- 
etration depth for an anisotropic superconductor with 
a sample consisting of a randomly oriented powder has 
never been attempted until now. 

NMR and electron spin resonance (ESR) have been 
used previously to investigate the anisotropy of MgB2. 
Two different components of the resonance signal have 
been identified in the superconductive state in a re- 
stricted range of magnetic fiel d 17 ' 18 ' 19 and the anisotropy 
of the upper critical field has been deduced. Addition- 
ally, an attempt was made to determine the temperature 
dependent penetration depth from the NMR linewidth 18 
assuming that MgB2 is isotropic, which is clearly not the 
case. 

Here we report n B NMR measurements in the tem- 
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FIG. 1: 1X B NMR spectra of the central transition in mag- 
netic fields of 1.97 T (left) and 3.15 T (right) obtained from a 
frequency sweep described in the text. The (±3/2 <-> ±1/2) 
satellites are about 340 kHz away from the central transition 
and, although they are not shown here, they are shown in 
Fig.6. 



perature range from 5 K to 40 K on a powder sample of 
MgB 2 at two magnetic fields, 1.97 T and 3.15 T. We find 
that on cooling the spectra acquire a broad asymmetric 
line below the superconductive transition temperature as 
shown in Fig.l. The shape of the broad line suggests the 
expected lineshape from an inhomogcneous field distri- 
bution from vortices in their solid state. However, as we 
will see, this interpretation is too simplistic. The relative 
weight of this broad line, compared to the narrow normal 
component, increases with decreasing temperature. We 
associate this with the temperature and angular depen- 
dence of the upper critical field. From this behavior we 
obtain the upper critical field anisotropy to be 5.4 at low 
temperature. We have also simulated the full spectrum 
at 8 K in a field of 1.97 T using this value for anisotropy 
and, by comparing with experiment, we have obtained 
the penetration depth A = 1, 152±50 A. Furthermore, we 
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FIG. 2: Schematic H C 2 diagram of MgE>2. The upper dotted 
line denotes 3.15 T, and the lower line corresponds to 1.97 T. 
The value for H C 2(0) is taken from Bud'ko and Canneld^S. 



find that the penetration depth is isotropic for T < 10 K 
even though the upper critical field and the coherence 
length are not. Our results support the theoretical claim 
that there are two different anisotropy parameters for 
upper critical field and penetration depth 3 -^. 

The polycrystalline MgB2 sample was prepared by 
solid state reaction techniques using a mixture of magne- 
sium and boron powders. The superconductive transition 
temperature was measured to be 39.5 K for the onset of 
diamagnetism in a magnetic field of 1.0 mT and 39 K 
for zero resistance. A sample of 0.2 gram randomly ori- 
ented MgB 2 powder was used in our experiments. NMR 
measurements were carried out in the temperature range 
between 5 K and 40 K in magnetic fields of 1.97 T and 
3.15 T in a superconductive magnet. Broad spectra were 
obtained by summing Fourier transforms of echo signals 
for a suite of different frequencies that cover the NMR 
spectrum. 

The spectra displayed in Fig. 1 are the central tran- 
sition (—1/2 <-> 1/2) of n B. At high temperature, the 
sample is metallic in the normal state and this spectrum 
consists of a single narrow and symmetric line. As the 
temperature is lowered, a broad and asymmetric line ap- 
pears. We associate this with the inhomogeneous field 
distribution from vortices in the superconductive state in 
addition to diamagnetic screening currents^. The weight 
of the broad line increases with decreasing temperature, 
while that of the narrow line decreases. The two lines 
coexist to a temperature of 5 K at 3.15 T, whereas only 
the broad line survives below 17 K at 1.97 T. This can 
be explained by anisotropy of the upper critical field in 
MgB 2 . 

Due to the temperature dependence and anisotropy of 
the two gap parameters in MgB^ 1 ^, its upper critical 



field has a temperature and angular dependence: 

H c2 (9, T) = Hg(T)/ y/l + tf B -l) cos* (1) 

where 9 is the angle between the applied magnetic field 
and the c— axis of a crystal. The temperature dependence 
of H^ 2 and is sketched in Fig. 2. For temperatures 
below T C (H), the upper critical field H C 2 will be equal to 
the applied magnetic field for crystals oriented at a cer- 
tain angle 9 cr {T). The crystals with 9 larger than 9 cr (T) 
have their H C 2 greater than the applied field, and are su- 
perconductive. Due to the random distribution of the 
orientation of the crystals, the superconductive fraction 
in the sample simply equals cos 9 cr (T) . As the tempera- 
ture decreases further, H^ 2 increases and, if it crosses the 
applied magnetic field, the whole sample becomes super- 
conductive. In Fig. 1, for H = 1.97 T, the narrow line 
disappears below 17 K leaving only the broad line. How- 
ever, in H = 3.15 T the c— axis upper critical field, H^ 2 , 
is always smaller than the applied field. Therefore, part 
of the sample remains in the normal state in this field 
and contributes to the narrow line in the spectrum, even 
at the lowest temperatures. 

In Fig. 1, the position of the narrow peak is al- 
most temperature independent and has a gaussian shape. 
Therefore, the contribution of crystals in the normal state 
can be deconvolved from the composite spectra. The ra- 
tio of the remaining area to the whole spectrum gives 
the superconductive fraction cos 9 cr (T), plotted in Fig. 
3. Furthermore, with H^O) taken to be 16 T from 
Bud'ko and Canfield 20 , and the upper critical field at 
9 cr at 5 K equal to 3.15 T, the external applied mag- 
netic field, the upper critical field anisotropy can 
be obtained from Eq. 1 and we find this to be 5.4 at 
low temperature. This value is consistent with previous 

Assuming to be temperature independent, H"^ at 
each temperature point can be obtained following Eq. 1. 
The temperature dependence for this analysis is plotted 
in Fig. 4 where it is compared with results for H"^ from 
other group o 19 ' 20 . The discrepancy grows with increas- 
ing temperature. However, it is now accepted that 
decreases with increasing temperature 3,7 . Therefore, in 
our derivation at high temperatures we have used a value 
for "fn that is too large which will produce a larger H"| 
and consequently an overestimate of the transition tem- 
perature in a given field. In contrast, Fig. 4 shows that 
the critical field curve deduced from our data and Eq. 1 
is too low. In fact, it extrapolates to a zero field transi- 
tion temperature around 30 K. The principal reason for 
this discrepancy is vortex dynamics. At high tempera- 
tures vortices are in a liquid stat o 16 i 23 and their dynam- 
ics on the NMR time scale average the local fields to zero 
at the U B nucleus. This transfers spectral weight from 
the broad line to the narrow line and reduces the appar- 
ent superconductive fraction obtained from NMR. At low 
temperatures, in the vortex solid state our analysis of the 
superconductive fraction is reliable and, as can be seen 
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FIG. 3: Temperature dependence of the fraction of supercon- 
ductive crystallites in the sample determined from the com- 
posite spectra plotted vs temperature for 1.97 T (uptriangle) 
and 3.15 T (downtriangle) . 



in Fig. 4, our results match H£|(T) below 10 K. 
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FIG. 4: H°2- The circles are results from F. Simon et alM 
(circle) and diamonds are from S. L. Bud'ko and Canfield^S. 
Analysis of our data using Eq. 1 is plotted assuming a con- 
stant 7h = 5.4 for 1.97 T (uptriangle) and 3.15 T (downtri- 
angle). We infer that our interpretation of the NMR signal 
as a vortex-broadened, inhomogeneous magnetic field distri- 
bution, is valid only below 10 K. 

We have found that the spin-lattice relaxation rate of 
the broad line is much slower than that of the narrow 
line. This agrees with a previous report^. Additionally, 
we have found that the rate increases smoothly with in- 
creasing frequency within the spectrum, rising in the high 
frequency tail of the central transition. Owing to the 
inherent inhomogeneity of the field distribution, which 
we will discuss later, it is not possible to deconvolutc 
spin-lattice relaxation signals to search for electronic ex- 



citations in different parts of the vortex structure as has 
been reported^ for YBCO. However, our spin-lattice re- 
laxation results serve as a guide to help us avoid selective 
saturation, particularly at low frequencies where the rate 
is small, allowing us to obtain a faithful representation 
of the spectrum. 

The absolute value of the penetration depth A is a key 
parameter for characterization of superconductivity and 
yet it is difficult to measure accurately. Using a tun- 
nel diode oscillator technique Fletcher et aL— found a 
penetration depth of MgB 2 between 800 and 1,200 A. 
Finnemore et al^ determined that A Q i, was 1, 400 A from 
transport measurements. Using ESR, Simon et al^ re- 
ported a value of the penetration depth between 1,100 
A and 1,400 A. From analysis of the second moment of 
the measured NMR linewidth, Lee et al. calculated the 
penetration depth to be 2, 100 A. But, as we mentioned 
earlier, the resonance methods cannot obtain a reliable 
measure of the penetration depth, if it is assumed in their 
interpretation that the superconductor is isotropic. Here 
we determine the penetration depth by comparison of our 
measured spectrum with a simulation of the local fields 
in the mixed state for an anisotropic random powder at 
8 K in a magnetic field of 1.97 T using the penetration 
depth as a variational parameter. 

The NMR spectrum is a local magnetic field map. 
At low temperature, the vortices are in the solid state 
and contribute to the associated field distribution of the 
NMR spectrum. The field distribution of the mixed state 
can be calculated by solving the Ginzburg-Landau (GL) 
equation. For this purpose we adopt Brandt's algorithm 1 
an iterative, quickly converging method. The solution 
gives the current and field distribution from the vortex 
lattice and the diamagnetic fields from screening currents 
in the superconducting state. The required inputs are the 
external field, coherence length £, and penetration depth 
A. We calculate the coherence length, £ a b = 108 A from 
the upper critical field 20 and we take its anisotropy from 
Eq. 1, 



t(9) = U/yfl + (j 2 H - l)cos 2 9 



(2) 



With these inputs, the field distribution for a crystal at a 
specific angle is generated by Brandt's algorithm 1 includ- 
ing the central transition and its quadrupolar satellites. 
We convolute the spectrum with a broadening function, 
exp— 2(H/S) 2 , which will also include the effect of the 
finite width of the NMR line in the normal state. In 
a powder sample, which we assume to be composed of 
single crystal ellipsoids of revolution, we must consider 
the shifts of magnetization owing to demagnetization ac- 
cording to the shape and orientation distribution of the 
grains 2 ^. For simplicity we characterize this distribu- 
tion by an average demagnetization factor, D. This as- 
sumption would be precise if the grain shape distribution 
is uncorrelated with the crystal structure. The demag- 
netization effect gives a relative shift of the magnetiza- 
tion which itself depends on the orientation of the grains 
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since the diamagnetic moment from screening currents 
is strongly angular dependent. Simulations of spectra at 
three different, but representative, angles are presented 
in Fig. 5. The spectrum for the whole sample is then 
obtained as the integral of 91 spectra with orientation 
uniformly distributed between and 7r/2, weighted by 
a factor sin 9 appropriate for a random distribution of 
grain orientations. 
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FIG. 5: Simulated spectra for crystals at different orienta- 
tions. The solid, dashed and dotted curves are the spectra of 
crystals with c— axis at n/2, n/4 and angles to the applied 
field, respectively. A demagnetization factor, D = 1/3, and a 
gaussian broadening parameter, 5 = 5.2 mT, were chosen for 
these spectra. 



of the first moment of the distribution, as can be seen 
in Fig. 5, must be correctly handled since it contributes 
significantly to the overall lineshape. Earlier work on 
other superconductors analyzing the field distribution in 
the mixed state has been directed at the moments of the 
distribution, so we have calculated the first three mo- 
ments for an anisotropic superconductor with randomly 
oriented grains, as a function of the penetration depth, 
restricted to the case of jh = 5.4 and 7a = 1. 
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FIG. 6: The spectrum at a temperature of 8 K and a mag- 
netic field of 1.97 T. The blue solid line is the experimental 
spectrum. The dotted line is the simulation described in the 
text. 



There are three variational parameters, A, D, and i5. 
We then carry out a \ 2 minimization of the difference 
between the simulated spectrum and the experimental 
one, taking their areas to be equal. The simulated spec- 
trum is shown in Fig. 6 together with the experimental 
spectrum. The numerical results provide an excellent 
representation of the complex measured spectrum with 
values for the variational parameters for the penetration 
depth A = 1, 152 ±50 A, average demagnetization factor, 
D = 0.31±0.01 and the gaussian broadening, S = 5.2 mT, 
which is larger than, but of the same order as, the normal 
state linewidth, 2 mT. The quoted accuracy is statisti- 
cal. This value for D is rather close to that anticipated 
for a spherical geometry, T) sp here — 1/3, and it is reason- 
able to expect this value for the average demagnetization 
factor for a large ensemble of grains. Earlier reports for 
the value of the penetration dept h 10 i 19 i 26 i 28 are similar 
to ours although our accuracy is higher. Our simulation 
and its comparison with experiment, as represented in 
Fig. 6, is the most precise such comparison obtained by 
resonance methods and it is the first time that such a 
simulation has been attempted for a strongly anisotropic 
supercondcutor. We emphasize that previous work has 
generally focused on moments of the measured spectrum, 
often restricting attention to the second moment. For 
an anisotropic superconductor the angular dependence 



The second moment of the magnetic field distribution 
of a spectrum from a vortex lattice can be related^ to 
its penetration depth A for low magnetic fields compared 
to H C 2 by the Pincus' formulaic where the second mo- 
ment varies as the inverse fourth power of the penetration 
depth, (B 2 ) = (0.0609</> ) 2 /A 4 . In the present case the 
simulated spectrum is the superposition of spectra with 
anisotropic coherence lengths and upper critical fields. 
Nonetheless, we find the 1st, 2nd and 3rd moments of 
the spectrum can be similarly related to inverse, even 
powers of the penetration depth in the following elegant 
way: 

(B) = -(l-D)-A 1 /\ 2 (3) 
(B 2 ) = S 2 /4 + A 2 /\ 4 (4) 

(B 3 )=A 3 /\ e , (5) 

where Ai, A2, A3 are numerical constants. The gaus- 
sian broadening factor is S and D is the demagne- 
tization factor. We find Ai = 1.415 x 10 4 T A 2 , 
A 2 = 2.621 x 10 7 T 2 A 4 , A 3 = 1.524 x 10 11 T 3 A 6 . 
However we caution that anisotropy and field dependence 
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FIG. 7: Simulation with different 7a plotted together with ex- 
perimental spectrum at a temperature of 8 K and a magnetic 
field of 1.97 T. The dashed lines are simulations. It is clear 
that the best comparison between experiment and simulation 
holds for 7a ~ 1. 

of the local field distributions mean that these numeri- 
cal constants hold only in a limited range which we have 
explored for MgB2 with 7# = 5.4, 7a = 1 and H = 1.97 
T. 

We have also investigated the effect of the penetration 
depth anisotropy 7a at low temperature. In actuality, 
the vortex structure for arbitrary angle 9 is found as the 
solution to the anisotropic GL equations—. However, as 
a reasonable approximation for an almost isotropic pene- 
tration depth, we introduce another variational parame- 
ter, 7a, and continue to adopt the solution of the isotropic 
GL equation for each crystallite. The anisotropy of H c i 



is the inverse of H C 2, therefore the penetration depth has 
the inverse angular dependence of the coherence length, 



\{6) = K^i + ^i-iy^e (6) 

With the same approach as before we generate the 
spectrum for the powder sample with 7a larger than one. 
The central transition is found to decrease with increase 
of 7a and the spectra become more asymmetric as Fig. 7 
illustrates. This suggests that 7a at low temperatures 
should be close to one. Magnetization measurements 
show that 7A is around 1.7 between 20 and 27 KM. SANS 
experiments on a powder MgB2 sample^ give an upper 
limit of 7a to be around 1.5 and essentially magnetic field 
independent. Our result is consistent with these values. 
However, SANS measurements on a single crysta&ii in- 
dicate that 7a is close to one at T = 2 K and at low filed, 
H < 0.5 T, and that it increases with external field reach- 
ing w 3.5 in a field of 0.8 T. The lower 7a in the powder 
sample is believed to be caused by a limiting crystallite 
size effect^. Further work will be required to elucidate 
this phenomenon. 

In conclusion, we measure the n B NMR spectra of 
a random powder sample of MgB2 in magnetic fields 
of 1.97 T and 3.15 T. The evolution of the spectra 
through the temperature range can be explained by the 
anisotropy of the upper critical field 7^, which is deter- 
mined to be 5.4 at low temperature. We find from our 
simulation that the penetration depth carries a differ- 
ent anisotropy from the upper critical field and that at 
low temperatures it is almost isotropic similar to that re- 
ported from SANSi£ for a powder sample. The value of 
the penetration depth that we have obtained for MgB2 is 
1, 152±50 A at 8 K in a magnetic field of 1.97 T. From our 
numerical studies we have found simple expressions for 
the penetration depth dependence of the moments of the 
field distribution in a random powder of an anisotropic 
superconductor. 
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